Hamiltonian dynamics of the two-dimensional lattice <^ 4 model 
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The Hamiltonian dynamics of the classical ip 4 model on a two-dimensional square lattice is 
investigated by means of numerical simulations. The macroscopic observables are computed as 
time averages. The results clearly reveal the presence of the continuous phase transition at a 
finite energy density and are consistent both qualitatively and quantitatively with the predictions 
of equilibrium statistical mechanics. The Hamiltonian microscopic dynamics also exhibits critical 
slowing down close to the transition. Moreover, the relationship between chaos and the phase 
transition is considered, and interpreted in the light of a geometrization of dynamics. 
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I. INTRODUCTION 



The study of the interplay between microscopic, deterministic dynamics and macroscopic statistical behaviour of 
large Hamiltonian systems is an old subject which dates back to Boltzmann. This subject has been mainly approached 
£j ■ having in mind the problem of the dynamical foundations of equilibrium statistical mechanics. In this framework 
the so-called ergodic problem is the central point: dynamics is studied in the perspective of proving ergodicity and 
mixing, thus giving a sound foundation to equilibrium statistical mechanics. Despite many efforts in ergodic theory, 
this goal remains distant, ergodicity and mixing having been proved only for abstract systems like the Sinai billiard, 
j-^. . However, a different approach is possible: instead of making a statistical assumption and then looking for its 
justification in the properties of the microscopic dynamics, one can consider dynamics from the very beginning. In 
practice, instead of considering a particular Hamiltonian system and trying to prove that its dynamics is mixing, one 
can observe the actual dynamical evolution of the system and measure the time averages of the dynamical observables 
of interest. In this way the statistical behaviour emerges directly from the dynamics, and one can wonder whether this 
behaviour is consistent with the predictions of statistical mechanics. In general, such an approach needs a tool that 
was not available at Boltzmann's times, i.e., a fast computer. In fact, this approach was pioneered by Fermi, Pasta, 
and Ulam fjj|, who performed, around 1950, the first numerical experiment on the relationship between dynamics and 
statistical mechanics, using one of the first computing machines, the MANIAC computer at Los Alamos. 

The conceptual point of view adopted in the present work is the Fermi, Pasta, and Ulam's one. We consider the 
dynamics from the very beginning, and we look at the statistical properties as emerging from the dynamics itself. The 
main goal of the present work is to show that a dynamical approach is worth and can provide some genuinely new 
understanding of phenomena that are usually treated in the framework of equilibrium statistical mechanics, like phase 
transitions. We argue that some of the information that is present in the dynamics of a Hamiltonian system, and that 
is thrown away at the beginning of the statistical-mechanical description, is relevant to the cooperative phenomena 
that show up in connection with the phase transition. This means that the relevance of dynamics in this context 
might go well beyond the foundational aspects related with the ergodic problem. 

In particular, by studying the dynamics of a paradigmatic system (belonging to the universality class of the two- 
dimensional Ising model), the </? 4 model on a two-dimensional lattice, we are able to show that the dynamical descrip- 
tion not only is consistent with the statistical-mechanical results, but also that a Hamiltonian description of dynamical 
phenomena like the critical slowing down is possible. Moreover, we suggest that a geometrization of Hamiltonian dy- 
namics based on simple tools of Riemannian geometry — originally introduced to describe chaotic dynamics |2|-|S|] 
- can provide a global description of the dynamical properties that are relevant to the statistical behaviour. We 
argue that most of these properties arc indeed the consequence of the actual geometric structure of the manifolds 
where the motion takes place. Hence geometry can be not only a useful tool in the theory of chaos, but hopefully 
can also provide the right language to bridge microscopic dynamics and macroscopic statistical behaviour which is 
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still lacking. In particular, phase transitions might be seen as the consequence of some major change in the geometric 
or even in the topologic structure of the "mechanical" manifolds [0] . Thus dynamics can bring in new concepts and 
tools in statistical physics. In our opinion these tools — here applied to a "standard" statistical-mechanical system — 
might prove useful both on conceptual and on practical grounds also dealing with "non-standard" topics in statistical 
physics, as the emerging field of phase transitions in finite systems (clusters, polymers, proteins) or long-studied but 
still unsolved problems like the dynamics and the statistical mechanics of glasses and more generally disordered or 
frustrated systems. 

Several other works have recently addressed the problem of the relevance of microscopic Hamiltonian dynamics to 
phase transitions, in particular in the framework of mean-field-like models p0[ and more generally as far as long- 
range couplings among particles are considered (ljj) . Moreover, there is now a renewed interest in microcanonical 
thermodynamics, both on general aspects (l^Jl^l and on phase transitions [ I^ , ^5| . 
The paper is organized as follows: in Section [n] we discuss some general 



aspects of the problem, then in Section 



[II the model and the relevant dynamical observables are introduced. Section IV is devoted to a discussion of the 
results of the dynamical simulations. In Section fvj we briefly recall the main points of the Riemannian theory of 
Hamiltonian dynamics and then we discuss the peculiar geometric properties of the two-dimensional ip 4 lattice model. 



Finally in Section VII we draw some conclusions. 



II. HAMILTONIAN DYNAMICS AND PHASE TRANSITIONS 

The dynamical aspects of equilibrium phase transitions are usually approached assuming from the outset the 
formalism of canonical equilibrium statistical mechanics. Dynamics is introduced only a posteriori, and usually 
making use of phenomenological, non-deterministic dynamics, e.g. Langevin dynamics |lq| . 

As already stated in the Introduction, our approach is different. Given a Hamiltonian system which exhibits a 
phase transition according to equilibrium statistical mechanics, we wonder what is its dynamical behaviour when it is 
studied as a natural dynamical system, i.e., associating a priori to the system a deterministic (Hamiltonian) dynamics. 
Accordingly, thermodynamic quantities will be obtained as time averages along the dynamical trajectories, whence 
- assuming ergodicity at least for the "usual" thermodynamic observables (see Ref. jl3| for recent results on the 
dependence of ergodic behaviour on the choice of the observable) — such quantities will be equal to microcanonical 
averages, whereas second-order equilibrium phase transitions are usually studied within the canonical ensemble. 

The two ensembles are equivalent only in the thermodynamic limit, thus the phenomenology observed in finite 
systems, as systems considered in numerical simulations necessarily are, might be different. To give only one example, 
let us consider the phenomenon of ergodicity breaking, i.e., the fact that ergodicity might no longer hold in the whole 
phase space but only in disjoint subsets of it. Such a phenomenon is indeed tightly related to phase transitions; in 
fact, when it happens, one may observe, as a consequence, a symmetry breaking, as in usual phase transitions. But 
ergodicity breaking is a more general concept than symmetry breaking, in fact one can recognize ergodicity breaking 
also as the origin of those phase transitions which do not correspond to the breaking of an evident symmetry of 
the Hamiltonian (for example in spin glasses) jl?]]. In the canonical ensemble, ergodicity can be broken only in the 
thermodynamic limit, while in the microcanonical ensemble, in principle, there might be ergodicity breaking also in 
finite systems. Ergodicity being a dynamical property, we think that a dynamical approach is particularly appropriate 
to study such a phenomenon. 

In the following we shall present in detail a dynamical study, pefrormed by means of numerical simulations, of 
the so-called </? 4 lattice model with Z2 symmetry. The first, preliminary, goal of such a study is to show that 
the dynamical phenomenology is consistent with the usual statistical description (the model exhibits a continuous 
phase transition which belongs to the universality class of the Ising model). The main results concern intrinsically 
dynamical properties, i.e., time correlation functions and Lyapunov exponents, and suggest interesting developments. 
In particular, two facts emerge: (i) the possibility of a Hamiltonian (i.e., ab initio) description of critical dynamics 
aspects like the critical slowing down, and (ii) a tight relationship between the local instability of the dynamics in 



phase space (characterized by the Lyapunov exponent and the related geometric observables discussed in Sec. VI) 
on the one side and the global phenomenon of the phase transition on the other side. It is worth noticing that a 
peculiar behaviour of the temperature-dependence of the Lyapunov exponent close to a phase transition (actually a 
Kosterlitz-Thouless transition) was observed for the first time by Butera and Caravati more than ten years ag o [ |18| . 
A more detailed discussion of the above-mentioned aspects as well as of the results can be found in Ref. p9[ . 
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III. MODEL AND DYNAMICAL OBSERVABLES 



Let us consider a discretized (lattice) version of the classical tp 4 field Hamiltonian, which reads as 

d 



i 



j 



\2 ^ 2 2,^4 

+3 - Vi) - 2 m + 41^ 



(1) 



According to equilibrium statistical mechanics, the system described by the Hamiltonian ([T]) has a critical point at 
finite temperature provided that d > 1. In the following we will restrict ourselves to the case d = 2; the case d = 3 
has been considered elsewhere, together with the cases of 0(n)-invariant ip 4 models po| . 



A. Dynamics and thermodynamic observables 

Our dynamical approach to the model defined by the Hamiltonian ([!]) is based on the direct solution of the equations 
of motion — Hamilton's equations — that read as 



tfi = Hi. 

d 



J £(<#+** + <Pi-» ~ 2( Pi) + ™?<pi - -yA<pf . (2) 



7T^ . \ j— / — |— j r i — j i — t i ! ■ ■ ■ —i 



The numerical integration of the 2N equations (|^) has been performed by means of a third-order bilateral symplectic 
algorithm pi] . The parameters have been chosen as follows: J = 1, m 2 = 2, A = 0.6. The average of any observable 
is defined as a time average, i.e. 



(/) = lim - f /Mr), n(r))dr ; (3) 

In practice, such an average is evaluated by means of a discrete sampling of /. 

As already discussed in Sec. [n], even if we cannot rigorously prove that the invariant ergodic measure associated with 
the Hamiltonian dynamics of our system is the microcanonical measure, nevertheless the microcanonical ensemble is 
the statistical ensemble which is naturally associated with Hamiltonian dynamics, for it is defined directly from the 
dynamics itself. Thus, in defining the dynamical observables which represent the thermodynamic properties of our 
model, we shall consider the microcanonical ensemble. 

The phase space density of the microcanonical measure can be written as [^2| 

ftnicro = -6(H(<p,n)-E) (4) 



where the normalization lo is given by 



J 5 (H((p, 7r) - E) dipd-n , (5) 



and dip dn is a shorthand for d<pi ■ ■ ■ dip^diri ■ ■ ■ d-K^. 

The entropy S is defined once the external parameter E (the energy of the system) is assigned. Different definitions 
of S can be given, all of which are equivalent in the thermodynamic limit. The two common definitions are (we set 
k B = l) 

S n (E)=logQ(E) , (6) 
S»(E)=logu>(E) , (7) 

where 

n(E) = f d (H(tp, 7r) - E) dip d-K . (8) 



The temperature is defined thermodynamically as 
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T dE ' 1 j 

in the above definition, S can denote either 5 or S u , whence we have two different definitions of temperature, T 
and T u . Again the two definitions are equivalent in the thermodynamic limit. Let us now limit ourselves to the case 
of natural dynamical systems, i.e., to Hamiltonian systems with N degrees of freedom whose Hamiltonian can be 
written as 

N 2 

H = y £g- + V(<p)=K + V; (10) 

i=l 

then the temperature can be expressed in terms of the kinetic energy as follows 0j : 

T° = - = !<if); (11) 

UJ 1 V 



- 1 I (K- 1 ) 



(12) 



It is worth noticing that for natural dynamical systems the definition (|llj) also coincides with the definition of temper- 
ature in the canonical ensemble. In what follows, unless explicitly stated otherwise, we shall adopt this definition of 
temperature and we shall drop the superscript f2, i.e., T = T n . A more general approach to microcanonical thermo- 
dynamics |H| has been recently introduced, which allows to find new expressions for thermodynamical observables. 
Anyhow, the definition (O) is the best suited for numerical simulations in the case of natural systems 

Another thermodynamic quantity to be computed in our simulations is the constant-volume specific heat, defined 

as 

j_ _ mm (13) 

C v ~ dE ■ [16 > 
From the definition and Eq. (Ill]) it follows that pjJ 

c v = ^-=[N-(N-2){K){K- 1 )]- 1 ; (14) 

in the thermodynamic limit N — > oo this expression for the specific heat reduces to the well-known Lebowitz-Percus- 
Verlet formula {SJ, 

which is almost universally used to compute the specific heat in molecular-dynamics simulations |26| . Notice that Eq. 
( |T5| ) is derived from an asymptotic expansion of the microcanonical fluctuations in terms of the canonical ones, thus 
it is valid only in the limit of large N. On the contrary, the formula (|l4| ) is exact at any value of N. Hence this is the 
correct expression to be used in a finite system. 

Finally, we will consider the order parameter, i.e. the "magnetization" 



(<P) = T?(Y.<Pi) ■ ( 16 ) 




IV. RESULTS OF THE DYNAMICAL SIMULATIONS 



We now turn to the problem of detecting the phase transition in our lattice ip 4 Hamiltonian system. According to 
the thermodynamic definition, we must look for a singularity of the thermodynamic observables as functions of the 
energy — or better as a function of the energy density e = E/N which remains finite as N — ► oo and facilitates the 
comparison of the results obtained at different lattice sizes. In particular, we look for a singularity in c v (e). Besides, 
being the transition associated with the spontaneous breaking of the Z2 symmetry, we will look for the appearance 
of a nonzero order parameter, i.e. for a nonvanishing value of the average magnetization (ip). 
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A. Binder cumulants 



In the canonical ensemble, a phase transition may occur only in the thermodynamic limit. As long as N is finite, all 
the thermodynamic quantities are regular functions of the temperature, and ergodicity and symmetry are not broken. 
Nevertheless, some marks of the transition clearly show up also in a finite system. The specific heat does not diverge, 
but exhibits a peak — whose height grows with the size of the system — at a temperature Xj?" (N) , while the order 
parameter is expected to be vanishing in the whole temperature range at any finite value of N. Neverthless this is 
true only in principle: in practice, e.g. in a canonical MonteCarlo simulation, where the simulation time is necessarily 
finite, the system is trapped in one of the two phases for a "time" which grows exponentially with N p7[, and one 
observes a fictitious symmetry breaking at a temperature Tf (N) . The latter temperature in general does not coincide 
with T C C "(A), even if 

lim T?(N) = lim Tjf(N) = T c °° . (17) 

In the microcanonical ensemble ergodicity breaking may occur also at finite N, hence one could expect a "true" 
critical energy to be defined also at finite N. No rigorous theoretical result is at disposal as far as this aspect is 
concerned. Nevertheless, it is reasonable to expect — and this is indeed what is observed — that the actual behaviour 
of the thermodynamic functions will be similar to the canonical case, at least as ./V is sufficiently large. In particular, 
we expect the specific heat to exhibit a peak at a critical energy density which is a function of N, and the order 
parameter to be nonvanishing below another critical energy, again depending on the size of the system. 

In the framework of the statistical theory of critical phenomena, finite-size scaling j27J] allows to estimate the critical 
properties of the infinite system from the values of the thermodyamic observables in finite samples of different sizes. 
In particular it is possible to locate the critical point by means of the so-called Binder cumulants p8[ . The Binder 
cumulant g is defined for our system as 

S = 1 "3W' < 18 > 

In the disordered phase the probability distribution of the order parameter will be nearly Gaussian with zero mean, 
hence g — 0. At variance, at zero temperature (or energy), when tpi = ip with no fluctuations, g — 2/3. At different 
sizes of the system, g{T) will decay from 2/3 to with different patterns. It is remarkable that the value of g at 
is independent of N, provided N is large enough for the scaling regime to set in, hence the critical point is given by 
the intersection of the different curves g{T) for different values of N. In principle, two different sizes are sufficient 
to locate a transition; in practice, owing to the unavoidable numerical errors which affect g, it will be necessary to 
consider three or more values of N. Moreover, the value of g at the critical point, usually referred to as g*, is a 
universal quantity; for a simple proof see e.g. Ref. [p7| . 

The theory behind the Binder cumulant method is totally internal to canonical statistical mechanics: to our 
knowledge, no extension of this theory to the microcanonical ensemble exists. Nevertheless we will adopt the pragmatic 
point of view of assuming its validity as a numerical tool also in our dynamical simulations, and our operative definition 
of the critical energy density will be the intersection point of the curves g(e) at different N . The consistency of 
the method will be checked a posteriori. In the following, unless explicitly stated otherwise, s c and T c will denote 
respectively and T c °°. 

The results for g(e) at different sizes for the two-dimensional lattice <p 4 model are shown in Fig. |l|. The crossing of 
the various curves at e c ~ 21.1 is evident. 



B. Thermodynamical observables 

1 . Temperature 

The temperature of the two-dimensional ip 4 system, numerically determined according to Eq. ([n]), is plotted in 
Fig. |2| Notice the change in the convexity of the function T(e) at e = e c . 

In Fig. H a comparison between the temperatures computed according to the two definitions (|ll|) and ( |l2] ) is shown. 
It is evident that already in a 10 x 10 lattice the two temperatures are practically identical. 

In Fig. |] the Binder cumulants are plotted vs. temperature T to locate the critical temperature T c . We see that 
this T c is consistent with T(e c ). 
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2. Specific heat 



The specific heat of the two-dimensional ip 4 system is plotted against energy density in Fig. ||. The asymptotic 
values of the specific heat in the limits e — * and e — * oo are exactly known. In fact at low energies the anharmonic 
terms in the Hamiltonian can be neglected, the system behaves as a collection af harmonic oscillators and c v « 1. 
In the high-energy limit the quadratic terms in the potential are negligible with respect to the quartic ones, whence 
Cy w 1/2 + 1/4 = 3/4. At intermediate energy densities, a neat peak shows up whose position is close to e c . The 
height of the peak grows with N. 

3. Order parameter 

Let us turn to the behaviour of the order parameter (ip). In principle, a nonzero value of (tp) is the characteristic 
signal of the breaking of the Z2 symmetry, which is in turn a consequence of ergodicity breaking. In practice, as 
long as a finite system is considered, the situation is more subtle. We have already observed that in a canonical 
ensemble, where the temperature is fixed, ergodicity can be broken only in the thermodynamic limit, hence the 
eventual appearence of a nonzero order parameter is a consequence of the necessarily finite observation time. In order 
to get reliable results the standard procedure is then to compute (\tp\) rather than (tp); in this way one has a quantity 
whose average is always nonnegative, and that coincides with the "true" order parameter as N — > 00. Obviously in 
the symmetric phase this quantity, at finite N, is nonzero: its amplitude will decrease as 1/y/N. This is precisely the 
behaviour of our numerical results reported in Fig. ^. 

In the microcanonical ensemble the situation is more complicated, in fact ergodicity breaking is no longer forbidden 
at finte N, hence a nonzero value of (<p) might be either a finite-time artifact or a "true" signal of symmetry breaking. 
Anyhow, the mere observation of the result is not sufficient to discriminate between these two alternatives. The order 
parameter (the "true" one, whose absolute value is taken only after the average for graphical reasons) is reported in 
Fig. fj]. We see that at each N there is a value of e below which the symmetry appears to be broken, and the value 
of (tp) seems to move almost abruptly from zero to a finite value. These "critical" energy densities are the closer to 
e c the larger N is. 

C. Critical behaviour 

The classical lattice tp 4 model, whose Hamiltonian is invariant under a discrete Z2 symmetry, belongs to the 
universality class of the Ising model. 

Our goal is not to obtain a precise measurement of the critical exponents, but only to check the consistency of our 
dynamical results with the statistical theory. The fact that the tp 4 theory belongs to the Ising universality class is a 
great advantage, because the critical exponents of the Ising model in two dimensions have been computed exactly: in 
particular, the order parameter critical exponent is /3 = 1/8 and the specific heat one is a = (the specific heat has a 
logarithmic singularity). We can compare the numerical outcomes of our simulations with the predicted Ising values 
in order to check whether the dynamically simulated critical behaviour is compatible with the predictions of statistical 
mechanics. In Figs. || and || the scaling behaviours of the order parameter and of the specific heat are compared to 
the exact Ising behaviours in two dimensions. The results are clearly very well compatible with the theory. This is a 
sign that the dynamical approach effectively reproduces the thermodynamical phase transition of the tp 4 model. 

V. DYNAMICAL PROPERTIES 

Up to now we have shown that the outcomes of the dynamical numerical simulations are perfectly consistent, both 
qualitatively and quantitatively, with the theoretical expectations regarding the phase transition of the tp 4 model. 
Though obtained through dynamics, all these results deal with equilibrium time averages: the time variable, even 
if not eliminated from the very beginning as in the statistical approach, has been integrated out in the averaging 
procedure. Now, we can also wonder whether there are intrinsically dynamical properties of our system that are 
relevant for the phase transition itself. Moreover, since from our point of view ergodicity breaking has its origin in 
the dynamics, we can try to understand what features are associated to a Hamiltonian ergodicity breaking. 
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A. Time correlation functions 



The first dynamical property that we are going to study is the dynamics of the order parameter tp(t). A qualitative 
understanding of what is going on is already provided by the time series ip(t) itself: some examples are reported in 
Fig. for various values of e, in the case of a 20 x 20 lattice. But a much more interesting information is contained 
in the time correlation functions of the order parameter, 



(<f (t)<p(t + T)) 



c v (r) = ^TJ;Z 11 ; (is) 



some of these functions are plotted in Fig. |ll|. We immediately note that close to the critical energy — and indeed 
very close to the "finite N" critical energy where ergodicity appears to be broken according to Fig. |?j — the shape of 
the correlation function changes rather sharply from an oscillatory pattern with a superimposed decay to a pattern 
indicating that the values of tp(t) are correlated over an extremely long period of time. Such a phenomenon is 
obviously reminiscent of critical slowing down. The latter is a feature of the localQ phenomenological dynamical 
evolutions, e.g., Monte Carlo dynamics, constructed in order to have the Boltzmann distribution as limiting 

invariant probability distribution, when the system described by 7i is close to a continuous phase transition. The 
usual explanation for the critical slowing down is that it is a consequence of the divergence of the static correlation 
length at the phase transition. However we are now considering a deterministic Hamiltonian microscopic dynamics 
and not a phenomenological MonteCarlo dynamics, hence our results clearly show that some kind of critical slowing 
down exists also in the microscopic natural dynamics. Something similar, i.e., the development of low-frequency 
collective oscillations, was observed in a planar mean-field Heisenberg model close to criticality fiofl . This result 
suggests that a Hamiltonian description of critical slowing down is possible and it might be very useful in dynamically 
understanding the phenomenon of ergodicity breaking (the study of a simple model which provides a first step towards 
this Hamiltonian approach to critical slowing down is presented in Ref. jl9| ) . 

In order to obtain a synthetic information from the time correlation functions, let us define a characteristic time r 
as follows: 

to 

C v {t)dt, (20) 



o 

where to is the time where C v has its first zero. Such a definition is to a large extent arbitrary, nevertheless it 
provides a relevant time scale either if the correlation function is oscillatory with typical frequency u> and with only a 
weak damping — as it happens at low energy — in which case r ~ u)~~ , or in the case of an exponentially decaying 
correlation with inverse time constant 7, in which case r » 7 _1 . The values of t computed with Hamiltonian dynamics 
are reported in Fig. [l2| The striking result is that the characteristic time is rapidly growing (notice the logarithmic 
vertical axis) as the system approaches the phase transition: the position of the peak is close to the "finite N" critical 
energy. 

B. Chaotic dynamics 

The lattice tp A model is a nonintegrable dynamical system. In the two limits e — ► and e — > 00, the system is 
integrable. The two integrable limits respectively represent a system of coupled harmonic oscillators and a system 
of independent quartic oscillators. The dynamics will be always chaotic in the whole energy range. Nevertheless, in 
analogy with other nonlinear oscillator systems, we expect that as the energy density is varied there exist different 
dynamical regimes characterized by different behaviours of the Lyapunov exponent A. In particular, the following 
questions naturally arise. Is there any peculiar behaviour of the Lyapunov exponent in correspondence with the phase 
transition? Is there any transition between different chaotic regimes in the tp lattice model, and, if yes, is there any 
relationship between these different dynamical regimes and the thermodynamic phases? 

There are not yet general and conclusive answers to these questions. Even if the study of a possible relation between 
chaos and phase transitions is a rather recent issue, which started with the already mentioned pioneering work by 
Butera and Caravati p8[, very different results have already appeared in the literature, ranging from the claim of the 



1 There are indeed "smart" dynamics which greatly reduce, or completely eliminate, the critical slowing down; the common 
feature of these dynamical rules is that of being highly nonlocal pOl] . 
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discovery of a "universal" divergence in A as the system approaches criticality in a class of models describing clusters 
of particles [|TJ , to the observation that the Lyapunov exponent attains its minimum in correspondence with the phase 
transition in Ising-like coupled map lattices |3^], to the apparent insensitivity to the liquid-solid phase transition of 
the Lyapunov spectra of hard-sphere and Lennard- Jones systems Jlll l ; for other recent results on first-order transitions 
see Ref. Q. 

More recently, some very interesting results have been obtained concerning mean-field-like models, in particular 
globally coupled rotators. Numerical results pi] , though of not easy interpretation due to strong finite-size effects, 
indicate that in these systems, that undergo a mean- field phase transition at a critical energy density e c , the Lyapunov 
exponent vanishes in the whole disordered phase, whereas it is positive in the low-energy (ordered) phase. Such a 
result has been theoretically confirmed in a very recent work [p5|. Since the latter work is based on the application 



of the theoretical tools described in Sec. VI, we will discuss there its results. 

Our simulation results are plotted in Fig. [l^. The first numerical evidence is that there is a strong dependence 
of A on N, what is peculiar of the presence of a phase transition. Moreover at large N a maximum of A(e) develops 
which eventually seems to move towards the critical energy density. Nevertheless no sharp, or singular, transition 
between different behaviours is found near e c , at variance with the three-dimensional case |20||; it is possible that such 
a sharp transition shows up as N — > oo, but no clear indications of this fact are provided by our results. Nevertheless 
the behaviour of A(e) in the region e < s c is very different from that of the thermodynamically disordered region, 
i.e, in the former A rapidly grows with e, while in the latter it has a quasi-flat shape (which is expected to become 
a decrease as e is large enough, because the high-energy limit is a integrable limit for the model). This behaviour is 
more clearly seen in Fig. 13 where a wider energy range is considered and logarithmic axes are used. This suggests 
that the phase transition has a dynamical counterpart in a passage between different chaotic regimes. However — at 
present — this statement cannot be formulated in a conclusive way, mainly because there is no clear and unique way to 
define and characterize a transition between different chaotic regimes (a recent improvement towards an unambiguous 
characterization of this kind of transitions can be found in Ref. p3| ) . 

Moreover, numerical experiments show that the detailed behaviour of the Lyapunov exponent close to the transition 
does not show "universal" features, i.e., it depends on the details of the Hamiltonian: on the contrary, in the low- 
energy range also for the tp model we have A oc e 2 that is the same behaviour observed in most, if not all, of the 
systems of coupled oscillators. Hence it is still unclear which feature of the e-dependence of the Lyapunov exponent 
has to be related with the phase transition. An exception is the already mentioned case of the mean-field rotators 
model p^l , where the order-disorder phase transition finds its counterpart in a chaos-order dynamical transition: this 
rather counterintuitive fact can be explained theoretically J35| and is likely to be a peculiarity of mean-field models. 
We shall come back to this issue in the following, in particular in Section VI when we will consider geometric properties 
which are strictly related with chaotic dynamics, and which, indeed, exhibit a much clearer behaviour near e c . 



VI. GEOMETRY OF DYNAMICS IN THE tp A MODEL 



Let us turn to the geometrization of the dynamics and its relations with the dynamical description of the phase 
transition. We will not enter the details recalling only the notations and the main results: all the details can be found 
in Ref. Q and in references quoted therein. 

The geometrical formulation of the dynamics of conservative systems was first used by Krylov in his studies on 
the dynamical foundations of statistical mechanics |36| and subsequently became a standard tool to study abstract 
systems in ergodic theory. Several new contributions to this subject appeared in the last years 

Let us briefly recall that the geometrization of the dynamics of ./V-degrees-of-freedom systems defined by a La- 
grangian C = T — V, in which the kinetic energy is quadratic in the velocities: 

' ^'u'l't ■ (21) 

stems from the fact that the natural motions are the extrema of the Hamiltonian action functional Sr = / Cdt, or 
of the Maupertuis' action Sm — 2 J T dt. In fact also the geodesies of Riemannian and pseudo-Riemannian manifolds 
are the extrema of a functional: the arc-length I = J ds, with ds 2 = g i jdq i dq' , hence a suitable choice of the metric 
tensor allows for the identification of the arc- length with either Sr or Sm, and of the geodesies with the natural 
motions of the dynamical system. Starting from Sm the "mechanical manifold" is the accessible configuration space 
endowed with the Jacobi metric 

(gj)ij = [E- V({q})} a l3 , (22) 
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where V(q) is the potential energy and E is the total energy. A description of the extrema of Hamilton's action Sh 
as geodesies of a "mechanical manifold" can be obtained using Eisenhart's metric |^7| on an enlarged configuration 
spacetime ({q° = t, q 1 , . . . , q N } plus one real coordinate q N+1 ), whose arc-length is 

ds 2 = -2V({q}){dq ) 2 + aij dq l dq } + 2dq°dq N+1 . (23) 

The manifold has a Lorentzian structure and the dynamical trajectories are those geodesies satisfying the condition 
ds 2 = Cdt 2 , where C is a positive constant. In the geometrical framework, the (in)stability of the trajectories is the 
(in)stability of the geodesies, and it is completely determined by the curvature properties of the underlying manifold 
according to the Jacobi equation [Q 

~d^ + R ^^ J IT- ' (24) 

whose solution J, usually called Jacobi or geodesic variation field, locally measures the distance between nearby 
geodesies; D/ds stands for the covariant derivative along a geodesic and R l jfem are the components of the Riemann 
curvature tensor. Using the Eisenhart metric (p3[) the relevant part of the Jacobi equation is 



^ + i?Vo^ = 0, i=l,...,N (25) 

where the only non- vanishing components of the curvature tensor are Roioj = d 2 V/dqidqj. Equation (pq) is the 
tangent dynamics equation which is commonly used to measure Lyapunov exponents in standard Hamiltonian systems. 
Having recognized its geometric origin, in Ref. || it has been devised a geometric reasoning to derive from Eq.(p5f) an 
effective scalar stability equation that independently of the knowledge of dynamical trajectories provides an average 
measure of their degree of instability. This is based on two main assumptions: (i) that the ambient manifold is almost 
isotropic, i.e. the components of the curvature tensor — that for an isotropic manifold (i.e. of constant curvature) 
are Rijkm = ko(gik9j m ~ 9im9jk), k = const — can be approximated by 

Rijkm ~ k(t)(gikgj m — gim9jk) (26) 

along a generic geodesic 7(^)5 (ii) that in the large N limit the "effective curvature" k(t) can be modeled by a gaussian 
and 5-correlated stochastic process. The mean fco and variance Ok of k(t) are given by the average and the r.m.s. 
fluctuation of the Ricci curvature kn — Kr/N along a geodesic: 

fc = (K R )/N ; (27a) 

al = {(K R -{K R )f)/N . (27b) 
The Ricci curvature along a geodesic is defined as 

1 „ dq l dq J 

where v 2 — ^-^jjr and Rij = R k ik j ^ s * ne Ricci tensor; in the case of Eisenhart metric it is 

K R ^AV = Y: -Jr- (29) 

i— 1 1 

The final result is the replacement of Eq.(|25|) with the aforementioned effective stability equation which is independent 
of the dynamics and is in the form of a stochastic oscillator equation |7j,|| 

^+fc(t)V> = 0, (30) 

where ip 2 ^ \J\ 2 - The exponential growth rate A of the solutions of Eq. (|30|), which is therefore an estimate of the 
largest Lyapunov exponent, can be computed exactly: 







where r = ■K^fk$/(2^Jko{kQ + <7fe) + 7rcr /c); i n the limit Ofc/fco <C 1 one finds 

A oc . (32) 

The latter result is the deep origin of the vanishing of the Lyapunov exponent in the whole disordered phase of a 
mean- field rotator model |34]; in fact, in that model, the curvature fluctuations can be analytically computed and 



turn out to be nonzero in the low-energy region but vanishing in the whole high-energy phase |35|, consistently with 
the mean- field character of the model. At the critical energy ou is discontinuous. The work reported in Ref. |3f| is 
particularly important because it is the first example in which the relationship between chaos and phase transitions 
is theoretically investigated within a framework which allows analytical calculations. The power of the geometric 
approach is evident in that case. 

It is natural to wonder whether the curvature fluctuations show any remarkable (singular) behaviour also in cor- 
respondence with non mean- field phase transitions. Numerical results suggesting a positive answer to this question 



have already been found in the cases of planar spin models || and three-dimensional 0(n) ip 4 models 20 1; in cor- 
respondence with a second order phase transition the curvature fluctuations show a cusp-like behaviour in all the 
cases considered up to now. These singular behaviours of the curvature fluctuations have been conjectured to be a 
consequence of a major change in the global geometry, if not in the topology, of the mechanical manifolds. 

In the following we are going to show that also in the case of the two-dimensional (p 4 lattice theory the above 
scenario is confirmed. Moreover, we will present some results concerning other geometric quantities, different from 
those defined in the framework of Eisenhart's metric; these results lend strong support to the topologic intepretation 
of the apparently singular behaviour of the curvature fluctuations. 



A. Curvature fluctuations with Eisenhart's metric 

In the case of the lattice (p 4 theory (|l|) the Ricci curvature per degree of freedom of the Eisenhart metric is 

fc* = ^=2dJ-m 2 + lf^. (33) 

i=l 

The time average k = (kn) of k^ is plotted vs. the energy density for the two-dimensional model in Fig. ^|. It is 
evident that k(e) changes its convexity close to e c , hence it shares this feature with the temperature. Anyhow in other 
models the shape of k{e) is completely different. 

What is more interesting, and, in the light of the above discussion, much more significant, is the behaviour of the 
fluctuations a\, reported in Fig. |16[ Actually in this figure the normalized fluctuation u^jk is reported. A cusp-like 
behaviour of the curvature fluctuations is evident in correspondence of the critical energy. 



B. Other geometric observables 

The same analysis of the time-averaged geometric quantities can be carried on also in the framework of the Jacobi 
metric — Eq. (J2^). In this case the ambient space is the accessible configuration space. We have studied the behaviour 
of the time average and of the fluctuations of the scalar curvature, whose corresponding dynamical observable is [p||lS|] 

U = A{E-Vf \- 2{E V)AV ~ {N ~ 6)1 W|2] ' (34) 

The results are reported in Figs, [I?] and |l8|. We see that also with the Jacobi metric the average curvature seems to 
have a smooth behaviour near e c , even if in the two phases the overall behaviours of the curvature seem different; on 
the contrary, the normalized fluctuation shows a sharp increase near criticality. Notice that these results have been 
obtained for small lattices only (up to 20 x 20) hence, in analogy with the other observables, we expect that in larger 
lattices this effect should be even more pronounced. 

Up to now geometry has been introduced through dynamics, by identifying the dynamical trajectories with the 
geodesies of suitable manifolds. 

However, other complementary approaches are also possible and interesting. Let us consider in particular the 
following one, recently introduced [ }40| in connection with the existence of different regimes in the dynamics of natural 
Hamiltonian systems. Given the dynamics, one can study the geometry of the trajectories as curves in the phase 
space, endowed with the Euclidean metric. In this way the trajectories are no longer geodesies of any manifold. 
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Nevertheless, being the Hamiltonian trajectories constrained on the constant-energy hypersurface the geometry 
of such curves carries informations on the "shape" of the invariant hypersurface. Such a relation can be made precise 
in terms of the geometry of He seen as a submanifold of R" Q. Let us define the curvature of the trajectory 
x{t) = ((fi(t), . . . , (fj\f(t), 7ri(t), . . . , 7Tjv(t)) in phase space as the generalization to n — 2N dimensions of the curvature 
of a curve in the two-dimensional Euclidean plane: 

(35) 

where r is the unit tangent vector to the curve x(t), i.e. r = £/|ir|, and s is the arc-length induced on x by the 
Euclidean metric of R™, hence ds/dt — |grad?i|. 

The time average and the fluctuations of the observable k can be defined as usual. In particular, the result for the 
normalized fluctuation (S 2 k)/(k) is reported in Fig. [l^ for the two-dimensional (p 4 model on a square lattice. We 
observe that the fluctuation of this curvature has a cusp-like behaviour in correspondence of the transition. Again we 
stress that this "quasi-singular" behaviour is obtained already with extremely small lattices. 

The phenomenology summarized above suggests that in correspondence of a phase transition the geometry of the 
manifolds underlying the dynamics undergoes a dramatic change. The precise nature of this change is still to be 
understood. Nevertheless, the fact that the singular behaviour of the curvature fluctations shows up using different 
geometric settings indicates that it is a consequence of some deeper property: the topological interpretation of these 
phenomena j^] is strongly supported. 



C. Geometric observables and Lyapunov exponents 



The geometric observables considered in Sec. |VI A| can be used to estimate the Lyapunov exponents. The result 
obtained applying Eq. (|3l]) to the two-dimensional ip 4 model are reported in Fig. |2(]. The agreement between theory 
and simulations is good qualitatively, but not from a quantitative point of view. A very good agreement is found only 
in an energy range of about two decades just above the transition. We do not have, at present, a deep explanation. 
A tentative one might involve the fact that in these models, especially close to the transition, the manifolds might be 
highly anisotropic (as witnessed by the growth of the curvature fluctuations) and thus the quasi-isotropy assumption 
might be no longer valid. Moreover, finite-size effects are expected to play a significant role. However, the non- 
satisfactory agreement between theory and simulation for the Lyapunov exponent is not a common feature of all 
models with phase transitions. In fact, in the case of the two- and three-dimensional XY models the geometric 
estimate of the Lyapunov exponent is in very good agreement with the simulations 0. Moreover, the analytical 
results obtained in Ref. [B5| compare well with simulations. In Ref. M] it has also been shown that in the case of ip 4 



models the estimate could be improved by adjusting the values of the timescale r in Eq. (|3l|). 



VII. CONCLUDING REMARKS 



The first result of the present work is that the statistical behaviour emerging from the microscopic Hamiltonian 
dynamics of a two-dimensional ip 4 model is perfectly consistent with the predictions of equilibrium statistical mechan- 
ics, and the agreement between the two descriptions includes critical behaviour. This means that all the tools coming 
from Hamiltonian mechanics can be used to investigate phase transitions. One of the aspects that can be considered 
in this perspective is certainly the possibility of an ab initio description of critical slowing down. 

A further result is that the study of intrinsically dynamical observables, as Lyapunov exponents, reveals an intriguing 
relationship between the local instability properties of the phase-space dynamics and the global phenomenon of 
the phase transition. In this perspective it is worth mentioning a recently proposed approach |4l| that unifies the 
description of Hamiltonian dynamics with the description of its stability through a path-integral formalism: this could 
be a way to relate stability properties with the phase transition from a general point of view. Here we approached this 
problem via a geometrization of the dynamics based on simple Riemannian tools. We found that the fluctuations of 
the curvatures of suitably defined manifolds associated with the Hamiltonian dynamics exhibit a "singular" behaviour 
at the transition. This fact is coherent with the results recently found for other models undergoing a phase transition, 
as the lattice ip 4 theory in three space dimensions, with both Z2 and O(n) symmetries pC|], and the classical XY 
Heisenberg model in two and three space dimensions 

A topological conjecture has been proposed to explain this behaviour || : the phase transition could be a consequence 
of a major topological change in the manifolds where the dynamical trajectories live. Such a conjecture receives further 
support from the results presented here, also because an apparently singular behaviour in the curvature fluctuations 
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is found also using different geometric settings with respect to those used in the previous works, which were all based 
on the Eisenhart metric. The problem of a precise characterization of these topological changes is still open and work 
is in progress in this direction (see Ref. jl9| for some preliminary results). 
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FIG. 2. Temperature T vs. energy density e. Symbols as in Fig. fil 
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FIG. 3. Comparison between the two different definitions of temperature T = T n — Eq. (|ll|), solid circles — and T — 
— Eq. (fL2|), squares — at three different lattice sizes: N = 10 2 (solid circles), N = 20 2 (circles), N = 30 2 (solid triangles). 
The solid line is the line T w = T n , the dotted vertical line marks the estimated value of T c (see Fig. GJ). 
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FIG. 4. Binder cumulants g (^) vs. temperature T. Symbols as in Fig. |l| The vertical dotted line marks the estimated 
value of Tc ~ 17.65. The inset shows a magnification of the transition region. 
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FIG. 5. Specific heat c v vs. energy density e. Symbols as in Fig. |lj 
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FIG. 7. Magnetization (ip) vs. energy density e. Symbols as in Fig. [I]. 
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FIG. 8. Scaling behaviour of the order parameter in the two-dimensional <p 4 model. Symbols as in Fig. [| only 20 2 , 30 2 and 
50 2 lattices are considered. The dotted line is the exact result for the Ising model, i.e., the power law \T — Tcj 1 ^. 
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FIG. 10. Temporal behaviour of the order parameter ip(i) for a 20 2 lattice at different values of e: from top to bottom, 
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FIG. 11. Time autocorrelation function of the order parameter <p(t) for the same lattice and at the same energies as in Fig. 
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FIG. 15. Average Ricci curvature k (^) with Eisenhart metric vs. energy density e for different sizes of the system. The 
symbols denote respectively N = 10 2 (solid circles), N = 20 2 (circles), N = 30 2 (solid triangles), iV = 50 2 (triangles). The 
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FIG. 17. Average scalar curvature (7V) (|3^) with Jacobi metric vs. energy density e for different sizes ol the system. The 
symbols denote respectively N = 10 2 (solid circles) and N = 20 2 (circles). The vertical dotted line marks the estimated value 
of e c . The inset shows a magnification of the transition region. 
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FIG. 20. Lattice cp 4 model. Synopsis of the geometric prediction of the Lyapunov exponent according to Eq. (^l|) (full 
squares) and of the numerical simulation results already plotted in Fig. [l4|. 
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